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_. Abstract 

C/^ , 

Q, ' We investigate the (reduced) Keller-Segel equations modeling chemotaxis of bio-organisms. We 

present a formal derivation and partial rigorous results of the blowup dynamics of solution of these 
equations describing the chemotactic aggregation of the organisms. Our results are confirmed by nu- 
merical simulations and the formula we derive coincides with the formula of Herrero and Velazquez for 
specially constructed solutions. 



> ■ 1 Introduction 

m 

^^ , In this paper we analyze the aggregation dynamics in the (reduced) Keller-Segel model of chemotaxis. 

, ' Chemotaxis is the directed movement of organisms in response to the concentration gradient of an external 

chemical signal and is common in biology. The chemical signals can come from external sources or they can 
be secreted by the organisms themselves. 

Chemotaxis is believed to underly many social activities of micro-organisms, e.g. social motility, fruiting 
body development, quorum sensing and biofilm formation. A classical example is the dynamics and the ag- 
gregation of Escherichia coli colonies under starvation conditions |17j . Another example is the Dictyostelium 
amoeba , where single cell bacterivores, when challenged by adverse conditions, form multicellular structures 
of ^ 10^ cells [inilSSI- Also, endothelial cells of humans react to vascular endothelial growth factor to form 
blood vessels through aggregation [22] . 

Consider organisms moving and interacting in a domain J7CM'',(i=l,2or3. Assuming that the 
organism population is large and the individuals are small relative to the domain fi, Keller and Scgcl derived 
a system of reaction-diffusion equations governing the organism density p : ft x M+ — > R+ and chemical 
concentration c : fi x R+ — ^ R+ . The equations are of the form 

dtp^DpAp-W-{f{p)Wc) 

dtc = DcAc + ap- (3c. ^ ^ 
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Here Dp, Dc, a, /3 are positive functions of x and t, p and c, and f{p) is a positive function modeling 
chemotaxis. Assuming a closed system, one is led to impose no-flux boundary conditions on p and c: 

d^p = and d,yC = on dQ, (2) 

where d,yg is the normal derivative of g. With these boundary conditions, the total number of organisms 
in a. is conserved. We refer the reader to [15l \T7\ [50l [58] for more information about chemotaxis and the 
Keller-Segel model. 

We presently concentrate on the case of positive chemotaxis, where the organisms secrete the chemical and 
move towards areas of higher chemical concentration. This leads to aggregation of organisms. Mathematically 
this is expressed as a blowup (or collapse) of solutions of ([T]). It was first suggested by Nanjundiah in [58] 
that the density, p, may become infinite and form a Dirac delta singularity. One refers to this process as 
(chemotactic) collapse. This is, arguably, the most interesting feature of the Keller-Segel equations. As 
argued below, the "collapsing" profile and contraction law have a universal (close to self-similar) form, 
independent of particulars of initial configurations and, to a certain degree, of the equations themselves, 
and can be associated with chemotactic aggregation. Though the equations are rather crude and unlikely to 
produce patterns one observes in nature or experiments, the collapse phenomenon could be useful in verifying 
assumptions about biological mechanisms|l| 

Phenomena of blowup and collapse in nonlinear evolution equations are hard to simulate numerically and 
the rigorous theory, or at least a careful analysis, is pertinent here. The recent years witnessed a tremendous 
progress in the development of such theories. We can now describe the shape of blowup profile and contraction 
law in Yang-Mills, cr-model, nonlinear Schrodinger and heat equations f [75 ] [73 l ISTl [5^ [70l [8l [56 l IST ] l29 ] ) [^. 
Yet, after 40 years of intensive research and important progress, we still cannot give a rigorous description 
of collapse in the Keller-Segel equations modeling chemotaxis. (See [I7l[HSIS71[Tni[IIl[ni[I31[IHl]for some 
recent works, [16], for a nice discussion of the subject, and [60l [45] [46] [42l [72] for reviews.) 

This is not to say that the Keller-Segel equations are harder than Yang-Mills, cr— model, or nonlinear 
Schrodinger equations, they are not, but neither are they less important. 

There are three common approximations made in the literature for system JT]) . Firstly, one assumes that 
the coefficients in ([T]) are constant and satisfy 



e := -^ « 1, 5 := — = 0(l) and f^ ■= ^ ^ 1- (3) 

The first of these conditions states that the chemical diffuses much faster than the organisms do. This is the 
case in practically all situations. As a result of this relation, one drops the dtc term in ([T]) (after rescaling 
time t — 7> t/Dp, this term becomes edtc). Secondly, one takes f{p) to be a linear function f{p) = Kp. Thirdly, 
the term /?c in ([T]) is neglected compared with ap, as one expects that it would not effect the blow-up process 
where p':^ 1 {\i is also small due to the last relation in ([3])). These approximations, after rescaling, lead to 
the system 

= Ac + p. 



^There are numerous refinements of the Keller-Segel equations, e.g. taking into account finite size of organisms f[T1[2l [55l l 
preventing complete collapse, which model chemotaxis more precisely. We believe the techniques we outline and develop here 
can be applied to these models as well. 

^ Numerical simulations for these equations failed until the compression rate was derived analytically, see [S] 1701 182) 
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with p and c satisfying the no-flux Neumann boundary conditions. 

Equations dH) in three dimensions also appear in the context of stellar collapse (see [40l [88l [24l [76] ) ; similar 
equations — the Smoluchowski or nonlinear Fokker-Planck equations — models non-Newtonian complex fluids 
(see [3TJ |53l [Ml EI] ■ This is the equation studied in this paper. 

We emphasize that in dropping the time derivative term of c, we have made the adiabatic approximation, 
in which the chemical is assumed to reach its steady state given by the second equation of ^ instantaneously. 

In this paper, we consider the collapse of radially symmetric solutions to the reduced Keller-Segel system 
([4]) on the plane R^ with a smooth, positive and integrable initial condition po and with the boundary 
conditions p, Vp, Vc — >■ as |a;| — ^ oo. To provide a right context for the discussion below, we mention that 
equation ([4]) has the following key properties: 

• It is invariant under the scaling transformations 

p{x,t) -^ -^P{j^' -yt) and c{x,t) -^ c{jx, —t). (5) 

• It has the static solution, 

• The total "mass" is conserved: /^^ p(x, t) dx — J^ p{x, 0) dx. 

We also mention that (Hj) (as well as (|T]) ) is a gradient flow, dtp = V • p\7£'{p), or dtp = — gradf (p), where 
£'{p) is the formal L^— gradient of £ and grad£'(p) is the formal gradient of £ in the space with metric 
{v,w)j :— — (v^J^^w) ,2- Here J :— V ■ pV < 0, whose inverse is unbounded operator, and £{p) is the 
"energy" functional given by 

^{P) = / (-^PA^V + Plnp - P) dx (7) 

(see Appendix [C] for more details) . We remark that the first term of £ can be thought of as the internal 
energy of the system and the remaining terms are the entropy. The solution ([B|) is a minimizer of £ under 
the constraint that J p = const . Note that J^2 Rdx = Stt, which is the source of Stt in ([8]). Under the scaling 
([5]), the total mass changes as 



/^p(^x,0)=A(-^)|p(x,0). 



Thus one does not expect collapse for d—\, and that collapse is possible for d > 2 with critical collapse for 
d = 2 and supercritical collapse for d > 2. (Equation Q in d = 2 is said to be L^— critical, etc.) 

Take po > 0. One has the following criteria for blowup of solutions of (|4|) ([59l[6]): If the dimension d = 2 
and the total mass satisfies 

M := / podx> Stt, (8) 

or, if d > 3 and ^f ^ po a; .^ sufficiently small (this means that po is concentrated at x = 0), then the solution 
to ([4]) blows up in finite time. 
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There is a fair amount of work done on equations ([T]) and (g]) and closely related equations. We give a very 
brief and incomplete review of it. Childress and Percuss [25] found that collapse for ([1]) with / linear does 
not occur when d = 1 and can occur when d > 3. For the two-dimensional case, they advanced arguments 
that collapse requires a threshold number of organisms. This threshold behaviour was confirmed by Jager 
and Luckhaus in [49] (see also [59l ES] [Ml ESI EH iMl IS^ ) ■ 

Herrero and Velazquez proved that there exist radial solutions of (|3]) for d = 2 with the threshold mass 
8|r2| collapsing to a Dirac delta singularity in finite time (see [37]). Also, unlike previous results, the authors 
give an explicit asymptotic expression of the developing singularity. They proved using matched asymptotics 
and a topological argument that for T > there exists a radial solution to (HJ) , which blows up at r = and 
t = T and is of the form 

1 f r<\{t) 

P(M)^^i?.,)(l + o(l)) + | ^^,-^,.(.-.1^) ^^,(^^^ 

as i — > T, where R\{r) :— R{r/X), R{r) is the stationary solution to (JU) (see (jl])) and 
\{t) = C{T - t)^e-75l '"(^-*)l^ I ln(T - t)\^\ i"(T-*)r^-i(i + o(l)). 

They also considered collapse of solutions to ([U with linear f{p) (see [38] and [39]). Obtaining similar results, 
they suggest that Jager and Luckhaus' adiabatic assumption does not affect the collapse mechanism. In the 
papers [54( i32j Lushnikov et al derived the log-log scaling as well as corrections beyond leading order log-log 
scaling. 

As noted in |37j . the asymptotics reproduced above are not of self-similiar type; that is, they are not of 
the form (T — t)^{r/(T — t)^) for some function $. In fact, as shown in [41], self-similiar blowup is not 
possible. Lastly, we mention that similar work has been done for the three dimensional case, where existence 
of collapsing shock waves has been shown. We refer the reader to [41] HO] [16] for these results. 

The above results are valid for radially symmetric domains and initial conditions. It was shown by 
numerous authors that the blowup threshold mentioned above decreases for the non-spherically symmetric 
situation. Moreover, Dirac delta singularities may develop on the boundary of the domain. We refer the 
reader to [5, 6, 33, 43l|47lET] for details. In 86|, Velazquez considers small radial and non radial perturbations 
of a collapsing solution and concludes, using formal matched asymptotics, that they are stable to these 
perturbations, leading only to small shifts in the blowup time and the blowup point. Existence of blowup or 
bounded solutions when f{p) is nonlinear was recently studied in ,48] . We also refer to j44) for a blowup result 
of a related Keller-Segel model. Lastly, we refer the reader to Horstmann [45l |46] for a more complete review 
of the literature including results on other models of chemotaxis and on the derivation of the Keller-Segel 
model as a continuous limit of biased random walks (see e.g. [67] Ell l80]). 

In spite of the considerable progress, the question of whether the mass collects in isolated points, forming 
Dirac delta distributions, remained unanswered. Moreover, these results give no information about the 
dynamics of blowup. These are the questions we address. 

Now, we describe the results of the present paper. Given a radially symmetric initial condition po{r) > 
sufficiently close to some R\„, for some Ap, and satisfying J po > J R, we show (formally, but with some 
rigorous supporting results) that the solution p{x, t) to (JH is of the form 

P(a^,i) = l^fiA(t)W(l + o(l)) 
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with A(i) — ^ as i — ?> T for some < T < oo. Thus, all the mass J pdx collapses to the single point a; = 
in finite time, or equivalently, the density p forms a Dirac delta singularity with weight Stt in finite time. 
Furthermore, we show that the compression scale. A, has the following explicit asymptotics 

X{t) = c(T - i)5e-75li"(^-*)l' I ln(r - i)|^(l + o(l)) (9) 

for some constant c. In Figure [T] we compare the blowup asymptotics ^ with direct numerical simulation 
of Q. 

We also give an estimate of the error term, p{x,t) — Rxu\{r), in the case when the nonlinear part in 
equation (1161) . given below, can be neglected. We believe that our results and our analysis can be made 
rigorous and can be extended to the full Keller-Segel system. 





Figure 1: The left pane shows the scaling parameter A(t) obtained by numerically computing the solution to 
(|4|) with the initial condition mo :— 4y^/(l + Sy + y'^)\y^r/\o with Ao = —6 — .01. The right pane plots the 

quantity J(A) := — (j^)^ against time, which according to ([9]) should be linear as the blowup time is 



approached. 



In^ 



We outline the approach used in this paper. In the case of radially symmetric solutions, the system Q, 
which consists of coupled parabolic and elliptic PDEs, is equivalent to a single PDF. Indeed, the change of 
the unknown, by passing from the density, p{x,t), to the normalized mass, 



m{r,t) ■■= — I p{x,t) dx, 

^'^ J\X\<T 



\x\<r 

of organisms contained in a ball of radius r, discovered by [JH [TB], maps two equations @ into a single 
equation 

dtm = A(°)to + r^^mdrm, (10) 

on (0, cx)) (with initial condition Too (r) '■— -^ J\ |< Po{x)dx). Here A" is the n-dimensional radial Laplacian, 

A" '■— r^("^^^9r''"^^9r ~ d^ + ^^^— !-9r. Thus (|4]) in the radially symmetric case is equivalent to (|10p and 
therefore we concentrate on the latter equation. 
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The properties of equation Q discussed above imply the fohowing key properties equation of PH)) 

• It is invariant under the scahng transformations m{r,t) —> m(j-r, y^^)- 

• It has the static solution (coming from the static solution R{r) = /-^ .^^.s^ of (U)), 

XW := YTr^- (11) 

• The total "mass" is conserved: 27rlim,._>.oo "^(j", i) — J p{x,t)dx — const. 

Note that the stationary solution has total mass 27rlim,.^oo x(^) = Stt, which, recall, is the sharp threshold 
between global existence and singularity development in solutions to (JH) (see ^). 

The properties above yield, as in the case of (|H), the manifold of static solutions A^o '■— {xi'''/^) I -^ > 0} 
and suggest a likely scenario of collapse: sliding along Aig in the direction of A —^ 0. To analyze the collapse, 
we pass to the reference frame collapsing with the solution, by introducing the adaptive blowup variables. 



r /■' 1 

m{r, t) = u{y, r), where y = t- and r == / —-— ds, 

A Jo A (s) 



where A : [0,r) -^ [0, oo), T > 0, is a positive differentiable function {compression or dilatation parameter), 
such that X{t) — > as i t T. The advantage of passing to blowup variables is that the function u is expected 
to have bounded derivatives and the blowup time is eliminated from consideration (it is mapped to oo). 
Writing (jlOp in blowup variables, we find the equation for the rescaled mass function 

drU = ^y^u + y^^uxdyU ~ aydyU, (12) 

where a := —XX. Now, the blowup problem for pO| is mapped into the problem of asymptotic dynamics of 
solitons for the equation (fT2|) , which was already studied in the pioneering works of [77l [78l [79l [TBI [13 IH31 

[HH [Ml [an [35]. 

The boundary conditions on u are 9"u(y, r) — > as y — > oo for a = 1, 2. As with the boundary conditions 
for Uni), these imply that mass is conserved: limj,_>.oo u(y, t) — limj,_>.oo u{y, 0). Equivalently, u, as a solution 
of ([T2l) . depends on a, which determines A, given A(0) = Aq, according to the formula 

X^it) = X^o-2 f ais)ds. (13) 

Jo 

Equation (fT2l) has the static solution (x(?/), a = 0). It is shown in [30] that the linearized operator on this 
solution has one negative eigenvalue —2a + t-^j + O (a ln~ i) (corresponding to the scaling mode — for a fixed 

parabolic scaling it is connected to possible variation of the blowup time) |j and one near zero eigenvalue, 
while the third eigenvalue, 2a + -p^ + 0(aln^ i), is positive, but vanishing as a -^ 0. (It also isolates the 

correct perturbation (adiabatic) parameter — -j — r;.) Hence we have to construct a one-parameter deformation 
of xiy) (besides the parameter A, or a). For technical reasons it is convenient to use a two-parameter family, 
Xbciy) 

XUV) ■■- ^, (14) 



''A similar analysis applies also in the subcritical case M < Stt where the solution converges to a self-similar one as t — ^ oo, 
which vanishes as f — > oo. In this case the operator Cah has strictly positive spectrum. 
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with 6 > 1 and both parameters b and c are close to 1, with an extra relation between the parameters a, b 
and c. The family Xbdy) gives approximate solutions to ([12]) (see (l43l) ) and forms the deformation (or almost 
center-unstable) manifold A4 := {xbd'''/^) I A > 0, p}. We expect that the solution to ([T^ approaches this 
manifold as r — > oo, and therefore we decompose the solution u{y, t) to (|12|) as the leading term, Xb{T)c{T){y) i 
and the fluctuation, (j){y,T), 

u{y,T) ^Xb{r)c{r){y)+(f>{y,T), (15) 

and require that the fluctuation 0(y, r) is orthogonal to the tangent space of A^ at Xb{T)c{T)iy) j {dpXp{T)i')j 0('i''')) 
0, where p :— (6, c). Note that this family evolves on a different spatial scale than <j){y,T) in (|15p . as it can 
rewritten as Xbcijj) ^ X>l,i{^) ^ Xbc- 

In parametrizing solutions as above, we split the dynamics of ^ into a finite-dimensional part describing 
motion over the manifold, M., and an infinite-dimensional fluctuation (the error between the solution and 
the manifold approximation) which is supposed to stay small. Substituting the decomposition (fT^ into the 
equation (|12p . we arrive at the equation 



^r<j)^~Cabc<i> + J'abc+^f{<t'), (16) 

where £ is a self-adjoint linear operator, J^ is a forcing term, and A/" is a quadratic nonlinearity. Due to the 
definition oi Ai, it turns out that its tangent space is very close to the subspace spanned by the negative 
and almost zero spectrum eigenfunctions (unstable modes) of the linearized operator, Cabc, and therefore (j) 
is (approximately) orthogonal to the latter subspace. 

The contraction law is obtained by using the orthogonality condition, (9bcXtic,0) = 0. The latter is 
equivalent to two conditions, 

dr{dbcXb(T)c(T){-),(t>{-.T)) =0 (17) 

and (9bcX6(T)c(T)(')' "/"('i '''))|t=o ~ Oi which lead, to leading order, to the differential equation 



2a^ 



(18) 



whose solutions, to leading order, are ^ (see Section!?]). 

We now describe the organization of this paper. In Section [2l solutions to @ are parametrized by the 
parameters (a, 6, c, (p) connected to u by (J12p . In SectionjH] we study the operator Cabc in (|16p and show that 
it has one negative eigenvalue and one simple eigenvalue near zero. We also give approximate eigenfunctions 
corresponding to these eigenvalues and prove that Cabc is positive on the space orthogonal to these quasi- 
eigenfunctions. In Section jj] we state the relationship between the blowup parameters a, b and c, whose 
proof is given in Appendix [D] and use it to obtain a dynamical equation for the blowup parameter a = —XdtX 
and derive the leading order behaviour of the scaling parameter A in terms of the original time variable. In 
Section (5] we derive the lower bounds on the operator Cabc- We use these bounds in in Section [6] in order to 
control the fluctuation (j) in the linearized equation, i.e. for ()16|) . with the nonlinearity M{(j)) omitted. 

In Appendix |X] we present the family of solutions to ((TT 



with mass 27r/i, where p, £ (2,4]. These solutions describe partial collapse with 2tt{4 — /i) units of mass 
concentrated at the origin. In the remainder of our work we will make no further use of these partially 
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collapsed solutions. In Appendix [B] we provide a proof of the orthogonal splitting theorem of Section [5] and 
in Appendix [C] we discuss te gradient structure of equations ^ and (JH) . 

In the following discussion, we use the notation f '^ g H there exists a positive constant C such that 
f < Cg holds. If the inequality |/| < C\g\ holds then we write / = 0(g). We also write / ^ g or / = o(g) 
if f{a)/g{a) — J- as a — )• and f ^ g ii the quotient converges to 1. 

Acknowledgements. The research of the second and fourth authors is partially supported by NSERC 
under Grant NA7901, and of the third author, by NSF under Grants DMS 0719895 and DMS 0807131. 

2 Parametrization of Solutions 

We parameterize solutions u\{y,T) of equation (|12p by the parameters a, b and c, and the fluctuation (j) 
according to 

u\{y,T) ^Xbc{y) + (l){y^T) (19) 

where A — A(t), b — h{T) and a = a{T). Substituting decomposition ([T^ into equation ([T^ gives that the 
fluctuation (j) satisfies 

a,0 = ~Cabc(i} + :Fabc + ^f{.4>), (20) 

where the linear operator, the forcing terms, and the nonlinear term are 

£,,, := -A(4) - ^^ - 1(6 - 1 - ^)dy + aydy, (21) 

(c + y ) y c + y^ 

,brC + bcr-2bca ,, 8(6 - 1) + 2ac - c^ 3260^(6 -1) 

Tabc ■■= -46, + 4 — ^ + 46c^ ' 2 2^3 ' 22 

c + y^ (c + y ) (c + y^)'^ 

AA(0) := y-V9,0. (23) 

Consider the weighted L^— space L^{R+,jabiy)y^dy), with the weight 



(c + y^)" 



7^ ^(y) = 7^^^' (24) 



and the corresponding inner product 

/"OO 

(/,5):-/ /(2/)5(y)7abc(y)y'dy. (25) 

The norm corresponding to this inner will be denoted by || • ||. The significance of this space is that, as we 
show below, the operator Cabc is self-adjoint on it. 

Remark. Another way to write Cabc is as 

Cabc = -A(°) - ^^^ - -^dy + aydy, (26) 

[c + y'^)'' c + y' 

and treat it as a self-adjoint operator on _L^(R+,7ab(y)y^dy), with weight j^j^J (y) = ,^1 'i)b , and corre- 
sponding inner product (/,.g) := J^ f{y)9{y) labc{y)y^^dy. 
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The decomposition ([T^ is not unique and as a result we have a single equation, (1^ . for four unknowns, 
a, 6, c and (p. Hence we supplement equation ((20|) with three additional equations. Two of the equations 
can be chosen as in [5H] to make the parameters a, 6, and c satisfy a chosen relation, say /(a, &, c) = 0. In 
addition, we have the relations 

(0,Cbc,>=O, 1-0,1, (27) 

in L^ {R'^ , -fabc{y)y^ dy) , for all times t > 0, where Cfeci are the tangent vectors to the manifold M := 
{Xbc{r/X) I A > 0, &, c}: 



j/^ 1 y^ 

o, ^--y^ucva/ . - — jT^, C&ci(y) := -rdbXbciy) = —- 2 



Cbco(?/) := -^^ydyXbciy) = ^„ , ,2^2 ' Cbci(t/) := -rdbXbciy) = 2 ' (2^) 



1 y2 

Cbc2{y) := --rrdcXbciy) = , , 2N2 ■ 
46 (c + y ) 

(The vectors Cbco{y) and Cbc2{y) are seen to be multiples of each other which confirms that one of the 
parameters is superfluous.) 

We proceed here differently and choose 

- (0, drCbci + [dr \^lab)C,bci) = - {Cabcfp, (bci) + {J^abc, Cbci) + (A/", (bet) , i = 0, 1. (29) 

As we will show the latter vectors are approximate eigenvectors of the operator Cabc having the negative 
and almost zero eigenvalues. In addition, we will choose a relation between the parameters a, 6, c. Eqns 
(pg]) imply that 

dr (0, Cbc^) = , (30) 

and therefore the inner products {(f>, Cbci) , * = 0, 1, are constant (one can think of this a constraint on a, b 
and c). The next proposition shows that ao,6o and co can be taken so that (Cbci,'/') |r=o = 0, and hence, 
by pop , we have (|27p . To be able to formulate a precise statement we introduce, for a fixed 6 > 0, open 
neighbourhoods of M , 

U,^{f: \\e-iy\f{y) - Xbc{y))\\oo < £ ,for some 1 < 6 < 2, ^ < c < 1} (31) 

and, for a fixed A > 0, 

Ue = {fir):fij)eU,}. (32) 

Proposition 2.1. Fix Aq > and < (5 <C 1. Then there is an e > and a unique C^ function g : lA^ ^>- 
((5, 1) X (1/2, 1) such that for mo G U^ we have the equation {C.bocoi^ 0o) = 0, or in detail, for i = 1,2, 

(too(Ao-) - XbocoXbocoi)\{ao,co)=g{mo) = 0. (33) 

For the proof of this proposition see Appendix [Bj 

Equations (|20p and (|29p form a system of coupled, partial and ordinary differential equations for the 
parameters (j), a, b, and c. We assume that this system has a unique local solution given initial conditions 
4)0, ao, bo and cq, the values of which are related to the initial value of m (recall that u\{y) = m{Xy)). 
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3 General Properties of the Operator Cabc 

Before proceeding we discuss general properties of the operator Catc mentioned above and used below, 

Proposition 3.1. The operator Cabc, defined on L'^{[0,oo),jabciy)y^dy) (with inner product (|25p J. is self- 
adjoint and has purely discrete spectrum. Moreover, we have the lower bound 

(<^,£afac0)>-[2a+ til±^]||0||2. (34) 

A/In- 

V ^ 

1 /2 

Proof One can check the self-adjointness of Labc directly or use the unitary map £_{y) -> 7^^^ (y)C(y)i from 
L'^ {[0, 00), jabc{y)y^dy) to L^ {[0, 00), y^dy) to map this operator into the operator 

Labc := llbc^^bclabc ^ (35) 

acting on L^{[0,oo),y^dy) with inner product (^,77) ;= J £,r]y^dy. The latter operator can be explicitly 
computed to be 



1„, „, 26(2 6-1 +ac) 46c 6 + 1 

-a^y^ - 2ab H !^^ '- - ——^ — --^ 

4 ^ c + 2/2 (c + y2 



iabc :- -AW + -a2y2 _ 2a6 + ^ \ , /^ '- - ]_ ' . (36) 



It is of Schrodinger type with the real continuous potential tending to 00 as j/ — >■ (X) as 0{y^). Hence, using 

standard arguments (see e.g. [36]), one can show that Labc is self-adjoint and its spectrum, and hence the 

spectrum of Cabc, is purely discrete. 

Now, we investigate the bottom of the spectrum of the operator Cabc- We begin with the operator 

Cf)bc '■— Cabc\a=(i- In what follows we use the convenient shorthand notation fx{r) — f{r/X), which we apply 

only for the subscript A. 

2 
Lemma 3.2. The function Cbc{y) '-^ 1 J' 2\2b is a zero mode of the operator C^bc '■= Cabc\a=o' 

\C-ry ) 

CobcCbc — 0. (37) 

The spectrum of £o6c starts with 0, which is a simple eigenvalue. 
Proof. Since xx is a static solution to (fTO| . differentiating the equation AV xa +''^^XA<9rXA = with respect 



C ■■- rdrx = j^-^ (38) 



to A at A = 1, gives that 
,2 
(1+72 

is a zero mode of the linearization of ()10p around x'- 

CoC = 0, (39) 

where Cq := — Av + r^^dr ■ X- (The vector Ca spans the tangent space of the manifold Mq :~ {xx | A > 0} 
at a point xa-) We deform this result as (|37p . Consequently, by the Perron- Frobenius theorem, the spectrum 
of Cobc starts with 0, which is a simple eigenvalue. D 
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The results above can be translated to the operator Li,c := Labc\a=o '■= loic ■^obdobJ i which is explicitly 
given by 



86c 46(1 - b)y ^ 



i>.-A'"-7r^ + iSSl^ m 



This is a deformation in b and c of the operator 

The operator Lbc, defined on the Hilbert space _L^([0, oo), y^dy), is self-adjoint with spectrum [0,oo). The 
bottom of the spectrum, 0, is a simple eigenvalue, 






LbcVbc = 0, with r]bc := ^loLCbc = -^--^Tb- (42) 



The tangent vectors Cbcoiv) and Cbciiv) (to the manifold Ai) are approximate eigenfunctions of the 
operator Cabc- Indeed, first observe that the functions Xbc{y) are approximate solution to p^ . Indeed, let 
$(w) be the map defined by the right hand side of (IT^ . $(u) :— Ay'u + y~^udyU — aydyU, then 

^, ,, , 8bca „, 4(6-l) + ac 326c2(6-l) 

c + y^ [c + y^y [c + y'y 

Now, differentiating $(x6c) with respect to c and 6 and using (|43p. we obtain 



,4( 6-l) + ac .^ c(6-l) 
c + y2 (c + j/^)^ 



/:a6caco(2/) = -2aCbco{y) + [4 ^ „ /J - 24f^^-^]aco(2/), 



^ . / ^ r 2ca 8c(26- 1),, , , 

CabcCbciiy) = [^-^ - ,\ ^j Kbciiy). 44 

c + y^ (c + y ) 

Though on the first sight Cbco and, especially, C,bci do not seem to be approximate eigenfunctions oi Cabc^ in 
fact they are. Indeed, assuming 6 — 1 = 0(a In -),c = 0(1), we obtain 

\\{Cabc + 2a6^o)Cbc^\\<lf'^^^'' J ^ J' (45) 

However, if one takes into account the normalizations 

||aco|| = -^ln^- + 0(^), ||C,,i(y)|| = ^ + 0(ln2i), (46) 

4V2 a 1^2 1 a a 

a 

then, for the normalized vectors we have 

. [- (b-l)+a j^Q 

|l(/:a6c + 2a5,o)^||<<^ V^ ' (47) 

llsfccill [a 1 = 1. 

The relation for i = implies ([M]) . D 
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4 Relation between Parameters a, b and c and Blowup Dynamics 

In this section, we state the relations between the parameters a, b and c, which is obtained by evaluating 
the equations in (|29p and is proven in Appendix |DJ Using these relations, we find the governing equation for 
air). 

Proposition 4.1. Let d := b — 1 and assume 

||</)||«(lni)-i, (48) 

" " a 

and, for simplicity, d < aln(a~^). Then 

Cr + Soicj), a, 6, c)ar = 4(^ - — y-) + O(j^) + 7eo((/), a, b, c), (49) 

dr + Si{4>,a,b,c)ar = -7-7^ + 0{-^ + a^d\n{-)) +TZi{(f>,a,b,c), (50) 

Ha) Ha) « 

where Si{(j),a,b,c) andTZi{(j},a,b,c), i — Q, 1, satisfy the estimates 

|5,(0,a,6,c)|<||0||(aln(a-i))'^\ (51) 

Remark 4.2. We see from (|52p that for the terms TZi{(j},a,b,c) in (j49p - (|50p to &e subleading we should 
have \\4>\\ < (aln(a~^)~i)i/2. 

As was mentioned above, this proposition is proven in Appendix [Dl Now, we choose a relation between 
a, c and d, so that the leading order term on the right hand side of the equation for c,-, (j49p . vanishes: 

d=-aln(a~^). (53) 



Proposition 4.3. Assume \\(f)\\ < . / . / °- 1 -, a?irf (1531) . T/ie^ i/ie function a{T) satisfies the differential 
equation 

which gives 

«W = ^(l + 0(r^)). (55) 

^T In r 

Proof. Plugging the relation ([53]) into ([50]) and remembering that d = & — 1, we obtain 

i(ln(a-i) - 1 + 2S2{(f>, a, 6, c))a, = -^^ + 0(^^) + 7e2(0, a, 6, c), (56) 

We see that to solve this equation for a-r, we need |52| ^ ln(l/a), which in view of ([5T|) with i = 2 requires 
that \\(j)\\ ^ ln(l/a). Due to the conditions of the proposition and estimate ((52)) with i = 2, the high order 
terms in equation (|54|) give a small correction upon integration and the leading part can be integrated exactly 
yielding ^. D 
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Remark 4.4. The above expression for a-r passes a consistency test: Or < and \ar\ <C a^. 

Proposition 4.5. For |T — i| <C 1, the scaling parameter X, with a = —XX satisfying (|55p . is asymptotic to 

A = fee- 3'"'^ (57) 

where r is related to t by 

fc(r-i) = -^e-5i"'^ (58) 

Proof. Using the definition a = —XX and tlie relation dtX — X^^drX we arrive at a = —X^^drX. Combining 

this with (|55|) we obtain the equation X^^drX — (l + 0(j-4— ))■ Solving this differential equation gives 

((57)) . Combining ([57]) with dtT — X^^ gives a differential equation for r(t) solving which in the leading order 
leads to (EHl. D 



Solving ([55]) for In r and substituting the result into ([57)1 gives ([S]). 

5 Lower Bound for the Operator jCabc 

In this section we investigate the linear operator Cabc, defined in ((2T|). The main result of this section is the 
following lower bound on the quadratic form (0, Cabc4>), <P -L Cbd, where, recall, the vectors (^bd are defined 
in 



Theorem 5.1. for |a| < 1, |fo-l|<l, \ar\ 't: a"^, {brl '^ a{l-b) and for any <p e H'^{[0,oo),jabciy)y^dy), 
(fiJ-Cbci, J = 0, 1 we have, for some absolute constant c > 0, 

{(f), Cabcff)) > ca ||(/'||^i . (59) 

Proof. Recall that the operator Cabc is unitarily equivalent to the operator Labc, 

Labc = llic^abclabJ^ , (60) 

acting on the space ^^([0, oo), y^dy) with the inner product (0, 77) :— J if)riy^dy. The latter operator has been 
explicitly computed in the proof of Proposition 13. II to be 

Labc^L,+Wiy)-2ab, (61) 

where 

(c + y2)2 c + 2/^ 4 ^ ^ 

and 

c + y' 

Since the lower bound of Labc is equal to the lower bound of Cabc, we estimate the former lower bound. 
We observe that, like Labc, the operator L^ is self-adjoint on i'^([0,oo), y^dy) and its spectrum is purely 
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discrete, provided a > 0. The latter property follows from the fact that the potential in (p^ goes to infinity 
as y — >■ c». Moreover, L* > 0. Indeed, write L* = Lofcc + \o?y'^, where 

(4) Ahc{h+l) 4b(b-l) 
Define 771(2/) :— ^ \ ) C(?/)i where C is defined in ([55|. so that £011^1 = 0- We compute 

An extension of relation (j63|) leads to the equation 

LobcVbc = 0, (64) 



where ?7bc is a deformation in b and c of 7/1 given by 
1 

(c-H/2) 



'^fc^ ■= 7^-:2\6- (65) 



Since Ty^c > we conclude, as in Lemma 13.21 that io6c > (with the zero being a resonance of Lbc)- This, 
together with L* = iofec + jCl^V^, implies that L* > 0. Next, we have 

Lemma 5.2. For <j) E H'^{[0, 00), y^dy), (pJ-d^bd, ^ = 0, 1, and a ^ 1, we have 

{L,^,^)>{4-0{^=))aml (66) 

'Ini 



where {(jj^rf) :— J (l)rjy^dy is the inner product in L'^ {[0, 00), y^dy) and \\ ■ ||* is the corresponding norm. 

Proof. To this end we will use the minimax principle for self-adjoint operators (see [74) ) . which states that 
the third eigenvalue, A3, of L* 

A3 = mf max , (67) 

where V is an arbitrary subspace of H^{M.*) and the estimate from 30^ of A3: 

A3 = 4a +-^(1 + 0(1)), (68) 

m - 

a 

for some constant C. Now, let 77 be the minimizer to {L^(j),(j)) over 
{(j)EH^'{\Q,^),y''dy)\dpLCbc^, * = 0,1, ||0|U = 1}. 

Since L^ is selfadjoint, 77 can be chosen to be real. Since the spectrum of L* is discrete, this minimizer exists. 
By the linear independence of C,bcii i = 0, 1 and orthogonality of 77 to Cbd, ^ = 0, 1, the three vectors rj, 
Cbci, 7 = 0, 1 span a three dimensional space. The minimax principle then asserts that 

A3 < max , (69) 

4>ew UK 
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where W :— spanJT;, Cf,cj, i = 0, 1}. Let (f>na, n = 0,1, he an appropriate orthonormal basis in spanjCbci, i = 
0,1}: 

(t>na ■= Vbci'na, (70) 

with 






(71) 



for sonie positive constants ci and C2- We write <p — 7177 + j24'0a + 730ia in the inner product (L*0, ^), 
where 

l7iP + l72p + l73p = l, (72) 

and use self-adjointncss of L* to obtain that 

(L*0, 0) =|7i|2(L^?7,?7) + \^2\'^{L^(f)Qa,(f>Qa) + |73 P (i*01a, 01a) 

+ 2Re(7i72)(77,i*(/>oa) + 2Re (7173) (77, L*</>iq) 

+ 2Re(7273*)(i*</'Oa,01a). (73) 

We compute the various matrix elements on the right hand side of equation (|73p . Using that 

A^^Vna = Ipna^y^Vbc + -qbc^yHna + 2(9j^f/'«a) (dyVbc) , 

we find 

i*'/'«a = ''PnaLobc'nbc + VbcHa-ipna ~ '^{dy'nbc) (dy'^Pna) ■ 



Using the facts that LobcVbc = 0, Haipoa = 2aV'oa and Haipia = 'iail^ia + {^C2aJ -^ - \/2cia)ij}Qa and 



computing 2(^dyrii,c) [dyipna) , we obtain that 

(L* - 2an)<j)na = Sn, 
with 

/„ /, N 2abc N , 
5o := - 2a 6 - 1 + ^^ '^oa , 
c + y 

and 

Si:^-2a{{b~l)-^^)(f>ia 
^ c + y^' 



- V2a{Ac2\\n-{h - 1) - i^Whi i + 2ci)</)oa . 
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Using this and the fact that the functions (j)ia are normalized, we estimate 

(i*0Oa,<^Oa)--^+O(^), (74) 

a 

{L.ci^oa, 0ia) = -^rr + o(^) (75) 

V21n- ln2 i 

fi a 

(L*0ia,<^ia) = -^(2 + ^) + O(-4T)- (76) 

ln2 i V2 In - 

a o. 

Let P^ be the orthogonal projection onto the orthogonal complement of the two vectors 0oa and (\)\a- 
We compute that 

||P^i*0.a||*=O(^— ), i = 0,l. (77) 

ln2 i 



Using the estimates, (l74|) - (|77l) . together with (j69l) and (fTS]). we find 

A3<max{|7iP(L,,7,,7)-^|72p + 0(^)}. 
7. 2 In^i'-' 

a 

Now, since L^ > 0, we know that [L^rj^rj) > 0. Then the above relation imphes 

A3<(i*ry,,?) + 0(^). (78) 

In 2 - 

a 

Using expression ([68l) for A3 in the last inequality, we obtain that 



(i,r,,77)>4(l-0(^^))a, 
/ini 



which gives the inequality (|66p . D 

Because of the decomposition (pT|) and since W{y) > and < fe — 1 < ,^ ^ , we arrive at 



/In -^ 

(</',ia6c<^)>^a||<^||^ (79) 

1/2 

or, by the unitary map (j) i—)- 7abc*^' 

{cf>,Catc<P)>laU\\\ (80) 

To pass from this bound to ([5^ . we decompose {(f>,Cabc4') = (1 — S){4>,^abc4') + 5{4'i^abc<P) and use (IM|) 
for the first term and Labc > — A^'*^ — C, for some C > 0, for the second one. Optimizing with respect to 5 
produces ((55|) . D 
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6 Analysis of Fluctuations 

In this section, neglecting the nonhnearity N{(j)), we find a bound < |1 ~ &| on the fluctuation </>. Given 
that we expect, from (|53)) . that 1 — b ^ fln^a' *^^^ ^^ sufficient to close the estimates. Neglecting the 
nonhnearity N((j)) in (I20p . we arrive at the linear equation 

dr(j) = -Cab(f> + J^ab- (81) 

More precisely, we have the following proposition. 

Proposition 6.1. Assume a, b and (j) solve (|8ip and (1291) (with M ~ Q), which is equivalent to (j27p . and 
are such that 1 — b — O(alni) and br — 0(a), and assume (|55p holds. Then, for t ^ 1, (j) satisfies the 
estimate 

II /ll / II //,^^ll / 2a , Ini ,1 

k^ ^0 r^ "+aln-. 
In i a 

a 

II II 2 

Proof. We use a Lyapunov argument with Lyapunov functional <j> i— >■ </> . The time derivative of this 
functional on solutions (p to ([5T|) is 

dr\\<j)f = -2 (0, /:,b,0) + 2 (J-,b„ 0) + (0, (5, ln7,b,)0) (82) 

We estimate right hand side of this relation. Let x(j/)? xiv) > be a smooth partition of unity, x^ + X^ = 1j 
s-t- x(2/) is a cutoff function that equals 1 on the set {ay"^ < k}, for some convenient large constant k > 0, 
and is supported on {ay^ < 2k}. We have 

Proposition 6.2. For any (p G H^{[0,oo),jabciy)y'^<iy),(l) _L C,bci,i = 0,1 we have, for some absolute con- 
stants ki,k2 > 0, 

2{(f),Cabc<t>) > aMh + /cia||0||^i + fc2(x0, a\'x0) ■ (83) 

Proof. Since is orthogonal to the vectors ^hd, i = 0, 1, and since a and b satisfy the conditions of Theorem 
IS.ll we have estimates ((59|) and (ISOp . which imply 



2{(j),Cabc^)>^aml2+kaml^. (84) 

Next, we estimate {(j),Cabc4') '^^ a different way. For the partition of unity defined after (1821) . we have the 
IMS formula (see e.g. [28]) 

Cabc = XCabcX + XCabcX " IVxl" - |Vx|'. (85) 

By ([34]) we have x^abcX ^ ^^X^- Using the inner product (^,77) := J £,iT^y and the notation := x'P we 
obtain 

{4> , Cabc<l>) = {(l>,llicLabcllic4>) , 

which, together with p6l) gives, for k large enough, 



i/2rl 2,2 o„^ 46c(6+l), 1/ 



/:a6c0) > (0,7:L[4«'r - 2«&~ (^;^2)2 ]7r'/') (86) 

>{^,(laY + ^Ka)^). (87) 
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Next, using that |Vx| and | Vx| are of the form ^/^X^ where x is supported between ay^ — k and ay^ — 2k, 
we compute |VxP + |Vxp ~ ^Xj which leads to 

('/',(|VxP + |VxP)<^><-||x</'lli.. (88) 

Using the IMS formula ([55)1 and the estimates above we find 

{<f>,Catc(f>) > -caWxc^Wh + ^(X0, (a'y' + '«a)X'^) " C^UWl^ , (89) 

for positive constants c, C. Now, write {(f),Cabc4') == (1 ^ 5){(l),Cabc4') + 5{cj),Cabc4')^ and use ([5^ for the 
first term on the right hand side and (|89)) for the second one, and choose 5 sufhciently small to arrive at 
dMl). D 

We substitute expression (|22|) for Fabc and observe that the orthogonality of to C^bd, « = 0, 1, implies 
1 A . / 1 



l,<h) =C( 7T,(t> ) =c 



^c + y^^ I \ (c + y^)^ ' 

to obtain 



Using Holder's inequality in the above equality implies 

(J-,b„0)=O((6-l) ||(j/>"-"0||). (90) 

Next, we estimate (</>, (9r ln7a;,c)0). Using that 9rln7a6c = — a^j/^/2 + 26,- ln(c + i/^) + 2&Ct-/(c + y^), we 
obtain 

(0, [dr \^labc)4>) = -^ar||#||i2 + 26,||(ln(c + y2))i/2</)||i. + 2&C,||^:^0||i. . (91) 

By ([50)1 and (I54p . we have Cr < 0, 6r < 0, and assuming c < 1, we have by (|49p and ((53)) . that Cr < 0. Hence 

(0,(9,ln7,bc)0> < -^ar||#|li2 . (92) 

Now, y^ < n/a on suppx, which implies (x0, (9r ln7abc)x</*) ^ "^S^llx^lli^- This, together with the relation 
(X'/'i (5t ln7Qf,c)x0) = (X0, (5r ln7a6c)X'/') + (x0, (^r ln7obc)X'/'), gives 

(0, {dr \njabc)cf>) < -^11X011^2 - ^arWymih ■ (93) 



Using the last estimate, together with (|82|) . ([SS]). (|90| . ([93|). (|53|) and (|55|) . we obtain, for some absolute 
constants fci, fc2, C > 0, 

a.||0||i. < -a||0||i. - fciallc/^lll,, - fc2(x0, (aV + '«a)x0) + C(aln i)||(y)-5+-</.||i. . 
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|2 



Using 9t-||</'|| = 2||(/)||9t-||(/)||, dropping the second and third terms (these terms can be used to control the 
nonhnearity) and dividing the resuhing inequahty by , we obtain 

drU\\<-'^U\\+Caln-. (94) 

2 a 

Now, integrating the last inequality gives that 

ll^ll <e-/o^'^('*)'^'*||0(O)||+ /" e--f>'''''^'^'{a\n-){cr)da. (95) 



/o 

We have computed that in the sense of asymptotic equivalence, a{T) ~ ^-^, as r — >■ oo (see equation 
(|55p). Consequently, as cr — >■ oo, we compute that 

r a{s) ds^^ Hra) In (^) = In (^) ^ '"^""\ 
Using 

e-/o"'^(^) = e--/'o'^°(")e" ■''>"('*) , 

noting that the first term on the right hand side is uniformly bounded, and the second term is ^ t^ts '"'^^ 
we obtain e~-'o "^"^^ — 0(r~T6 '"'^). Using now 

a 

we obtain that the term involvmg the mitial eondition m (j95p is bounded as 

e-/>(^)^^||0(O)||<||0(O)||(-^)'"^. (96) 

in - 

a 

To bound the integral term in ()95p we begin by splitting the domain of integration into [0, ar] and [ar, r] 
for some < a < 1 to be chosen later: 



1 /"T 1 



Since (aln-!-)(o') is decreasing and cr i— > e J" ''('*) '^^ is increasing and both are positive, we can bound these 
terms from above by 



i:.r<'^)ds [ (a\n-){a)da+{a\n-){aT) j e' ^^'^'^ '^'^ da 



Since e Jot "('')'''' ^ CaT^^"^"', the first term is bounded from above by Ca(aln i)(0)r^+2 '"". Taking a 
such that ln(Q;)/2 < —2 the first term is < r~^, and the second is bounded by (a In i)(aT) ~ Cln^ t/t. So 
we find 

e-J'-"'^(^)''^(alni)(a)dc7< (alni)(T). (97) 

Using bounds (IMl) and (^7]) in (^5]) completes the proof. D 
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A Complete Set of Static Solutions for the Radial rKS 

The static solutions of equation (TTOl) satisfy the second order differential equation 



d^,X+-{X'l)drX = (98) 

r 

and hence form a two dimensional manifold. We prove 

Proposition A.l. Equation (|98p has the one-parameter family of static solutions 

(and therefore the two-parameter family Xa ('") •= X {^/^) o,-'' well). The mass at infinity of Xx ^s M- 

Remark A. 2. If fi < A, then the mass at the origin is non-zero, i.e. blowup has already occurred. If fj, > A, 
then the mass at the origin is negative and hence the static solution is not physical. 

Proof. We use the transformation 

Vix) = r 

X 

in (j98p under the assumption that the right hand side is indeed a function of x alone. Using this transfor- 
mation, equation (|98l) becomes 

Xdxip + "0 = 2 - X. 
Integrating this equation gives that 

1 n 1 

^ 2'^^2x 

and hence, upon substituting this into the definition of ip and integrating over r, we obtain the general 
solution 

X 



l+(i)^ 



where r± — 2 ± \JA + v are the roots of x^ — 4x — J^ = 0. The total mass at infinity of these solutions is 
r+ and hence it is natural to define a new parameter /i = r+ G [2,oo). The static solution in terms of the 
parameters A and v are 

,_ (irv+4-/i 

The constant A is positive since it is the exponential of the constant obtained in the last integration. D 



KS July 31, 2012 21 

The tangent space of the manifold A/a.^ is spanned by the functions 

and 



_ y^^t"-^^ -l-2y^'-^\ny 

t^^^'t^ - (1 + j/M-2)2 ' 



where y = j. 

We again restrict to the situation of k = 1 and n = 4. After the gauge transform the zero modes Ca u 
and Ca u transform to 

''' '~ A2 h2 ''A.m ^"^^ '^A,M ■- ;^2 2 ''A.M' 

neither of which are in L^ (r'^ dr) and hence are generaUzed eigenfunctions of C (without the A term. By the 
ODE theory the above functions are the only linearly independent solutions to the equation Cq(J) — 0- The 
Perron-Frobenius theory shows that is the lowest point of the spectrum of C. 

B Proof of Proposition 12.11 

Both existence and uniqueness follow from a standard implicit function theorem argument. Fix < (5 <C 1 
and let Z := e^y L°°{[0,oo)). Recall that 6=1 + l/2alog(l/a), and define the vector-valued function 

G(/, a, c) := ( (/ - Xbc, Cbc^) , * = 0, l) (99) 



OO 



= (^ y (/(y) - Xbciy)) (c + y^'-'+^ye-^y' dy, z = 0, l) . (100) 

This function maps Z x M_|_ x R+ into M^. It is a C^ function and G{xbc, a, c) = 0. Moreover, the derivative 
of G with respect to (a, c) at / = Xbc is 

A:=( j;"" 1°^), (101) 

i la i Ic 

where 

r.a ■■=-fj y^{c + y^r-^+^e-^^y^Ay (102) 

r.c:=-^y^ 2/3(c + yY^-4+»e-t^^d2/. (103) 

Compute that the determinant of A satisfies 
|detA| = -l^(ln2i+0(l)), 
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as a -> 0, and so | det ^| > C > 0, for some constant C, for (a, c) G (0, S) x (1, 2) for S small enough. Thus 
by the implicit function theorem, for any a* G (0,5) and c* G (1;2) there exist open sets Ua,c, C Z and 
Ki.c, C (0, 5) X (1, 2) containing x&.c and (a*, c*), respectively, and a unique function ga,c, '■ Ua,c, -^ ya,c, ■ 
To determine the size of the neighbourhoods Ua,c, we look more closely into a proof of the implicit function 
theorem. Write /i = (a, c) and expand 

G{f,^l)^G{f,^l.) + ^^G{f,^l.){^l-^l.) + R}{^J), (104) 

where Rf{lj) = 0(|/i — ^u^p) uniformly in / e i?c(Xf).c. ) and {a^,,c^) £ {5/2,5) x (1,2), for any fixed constant 
C. By continuity and the above computations, there is e > such that dct 9^G(/, /i*) is bounded away from 
zero uniformly for / g Be{xb,c,) and {a*,c*) S {5/2,5) x (1,2). From ()104p we find a fixed-point equation 
for ^ — ^^, 

fj.- H^ = ^f{^i - Ai*), 

where 

1>/(Ai) = -{^Mf,^^*))~\G{f,^i.) + Rf{fi)) . 

The above observations imply that there is an ei > such that $/ is a contraction on Be^{tx^) for any 
/ G Be{xb,c,) =: Ua,c,- Taking the union of Uac over a G (5,1) and c e (1/2,1) gives the open set Ue- 
Patching together the functions Qac gives g. D 

C Gradient Formulation 

The Keller-Segel models ^ and @ are gradient systems. We begin by formulating a normalized version of 

9iP - Ap - V • (/(p)Vc) ^^^^^ 

e9fC = Ac + p — 7c, 

as a gradient system. This system is obtained from ([1]) by setting unimportant constants to 1. 
Define the energy (or Lyapunov) functional 

£f{p,c):^ j h^Vc\^-pc+^c^ + G{p)dx, (106) 

where G{p) :— f g{s) ds and g{p) :— f jj-p- ds. The L^-gradient of Sf{p, c) is 

grad,.£,(,,c)H_rc-p + ,c)' 
and hence, if we define U = {p,c), then (jlOSp can be written in the form dtU — I£'f{U), where 



V • /(p)V 




/ = l „ 1 



KS July 31, 2012 23 

The operator / is non-positive and may be degenerate, however, assuming it is invertible, the operator / 
defines the metric {v,w)j := — (w,/~^w),2 ,3. In this metric, gTad£{U) = —I£'{U) and hence 

dtU = -grs.dEf{U). 

This shows that (|105p has the structure of a gradient system. A consequence of this is that the energy 
decreases on solutions of the KS system. Indeed, if / > 0, then 



9tf/(p,c) = -||/(p)W(c-5(p))f --||Ac + p-c 



111. , l|2 

£ ' 



The gradient formulation for ^ is similar to the one for (|105|) . Instead of ()106|) . one uses the energy 
(III). The latter is obtained from ()106p by dropping the quadratic term ^c^, replacing c with — A~^/3 in the 
remaining terms and using that /(/?) = p. The formal Gateaux derivative of £ is dpE(p)(l) — J(A^^ p +\ap)4), 
and therefore the gradient in the metric (u, w) j := ~ (u, J~^wj ,21 where J := V • pV < 0, is 

grad£(p) = -V • pV(A"^ p + \np) = -V • pVA~^ p - Ap, 

which is the negative of the r.h.s. of the first equation in (j4]) with c — — A^^p. Hence the equation ([4]) can 
be written as dtp — — gradf(p) in the space with metric {v,w)j :— — (y, J^^vo) ^2- Again, the energy £ 
decreases on solutions of (|3]): 

dt£ = {E'JE') = -\\p'^VE'\\^. (107) 

This can be thought of as an entropy monotonicity formula. 

The stationary solutions of ^ are critical points of the energy functional £", given in ([7]), under the 
constraint that J p = const.. Thus, they satisfy S'{p) — C, where C is a constant. Explicitly S'{p) = C 
reads 

log(p) + ^P = C ^ Alog(p) + p = ^ AM + e" = 0, (108) 

where u = log(p). Solutions to (|108l) can be written in the form of 'Gibbs states' p = M-f—^ (see [33]), with 
the concentration c considered as a negative potential (remember that Ac = —p). In two dimensions, this 
equation has the solution R — -rrrr-Tryi- This solution is a minimizer of S under the constraint that J p = Stt. 

D Proof of Proposition 14.11 



In this appendix, we prove Proposition |4T] relating the parameters a,b and c, by evaluating the equations 
in dMl)- 

Proof of Proposition \4-. 1] Let Ri{(})) :— {Cabc4'^(,bci) — {J^,Cbci) ■ Here and in what follows i = 0, 1. The 
equations 1^^ can be rewritten as 



{J^abc, Cbci) + (0, drCbci + (dr l'ajab)Cbci) = Ri{(f>)- (109) 



c + y 
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We begin with evaluating {Tabc, Cbd) to leading order. To this end, we begin with the elementary computation 

{l,Cbc^)^2"-^a-"-^+0{^\n'-), (110) 

a' a 

^, Cbc^) = T-^a-' \n^-' - + 0{ln^' -), (111) 

I a a 

, \,^ XbcS = 2-5c'-i In^ i + 0(ai-» In^-' -), (112) 

yc + y^Y I a a 

( ^\,3 ,acA=2'-gc'-^ + 0(alnl). (113) 

(c + y^)-^ / a 

These estimates are proven at the end of this appendix. Using these estimates in ^2f2\ . we arrive at 

(^a&cCfcc.) = -6^2'-2a-»-i(i _ c2-ialni-' -) + cAhT"^ oT'' h^^" - - 6c'2"3ln' -) 

^ a a a 

- V'-Hca}'" In^-* - + (6cV + ah^"^^^-'^) W ^ - 6cM2*-i 

a a 

+ 0{bra-' In^ -) + 0{{ar + a) In^* -) + 0{{d + a + Cr)a^-' h?-' -) + 0{da In -) . (114) 

a a a a 

Next, we compute the term {(j),drC,bci + (c^t ln7abc)C6ci)- Differentiating Cf,ci and 97-ln7a& with respect to r, 
we obtain 

a.a« + {dr \nlabc)Cbc^ - [26. ln(c + y^) + (26 - 2i-*)^^ - ^y^]^^^^^ (^15) 

■- c + y^ 2 -■ 

Using in the case i = 1 that is orthogonal to CbcOj we find 

((/), 5rCbci + (9r ln7ab)C6cj) = ^rS'il (</>)+ Crfe (0) - 0x5*43 (0), (116) 

where 

5,;i(0):=2(0,ln(c + y")Cbc,:) 

Sa{(f) := 2^,o(6 - 1) U. -^ibc^ 
\ c + y^ 

Sa{(l)) := -j^{(t>,y'^Cbci) ■ 
Collecting pl4)) and (fTTB)) . we have, for i = 0, 

(5o3(0) + 0(l))ar + [^ - ^ ln(a"') - 5oi(<^) + 0{ln^ -)]br 
4a 8 a 

- [^(In(a-i) - 1) + S02W + 0{a\n -)]cr - -^[iacln(a-i) - d - ac] - i?o(0) ■ (117) 
and, for i = 1, 

(5i3(0)a + 0(ln2 i))a, + [-!-£- a5ii(0) + 0{a-' In" -)]br - [^ - ia6cln(a-i) + a5i2(</)) + 0(l)]]c. 
a 2a 4 a 4 8 



ll-dHa-^)-^ 



-abc[- - dln{a-^) ln{a-^) + d] - ai?i(<?!)) . (118) 
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We manipulate equations (|117p and (|118p and solve tlieni for br and c^ to obtain 

- foar + gcr ^ -bcln{a^'^) - -— + vo + ro, -fiar+gbr^ — + Vi+ri, (119) 

8 4a 8 

where 

1 c 1 c 

/o := (^ T - a5'ii)5o3 + (t + z"-Mo'^^) - aS'oi)S'i3, 

6 In(a-i) lb be 1,2,1 -1^2 1, , I -l^ .2 .c ^ c ^ 

16 a 4 2a 8 16 8 a 

- -ln(a"-^)(6cS'oi + 6511 - cS'12) +a(-65ii - S'o2S'ii + 5'oiS'i2), 

o 

vo := -6c[2dln(a"i) + 3 + d + caln(a"^)((iln(a"i)- --d)-ca- -a^c^\n{a^^f\, 
8 2 4 

ro := (- Tc)i?o - -i?i + caln(a"^)(6dS'oi - -i?i) 

+ a(-6cS'oi + -6^511 + i?oS'ii + S'oii?i) + -a^6c(5ii - - ln(a"^)(cS'oi + S'n)), 
/i :=a(-ln(a"^) - - + So2)Si^ ~ '-I ~ V"l'^(""^) + a'5'i2)S'o3, 

1 3 11 

Vi := 6 c[daln (a^ ) (ialn(a^ ) — a( d) H — ca ln(a^ )], 

8 2 2 4 

ri := --i?o + Kz^Ro + -dcSo2 + -i?i)aln(a"^), 
4 o o 

+ a(— 7; 1- bcdSQ2 H ;^ ^0512 - -bRi + R1S02), 

Z Z o 

+ lbca^S,2 + \cS02Ha'') + i5i2ln(a-i)). 

Next, we derive estimates on Sii{(f>), Si2{(f>), -Ro('/') and i?i (</>). Using the Cauchy-Schwarz inequality 
and simple modifications of the estimates P^ , we arrive at the estimates 

\S^im < ||0||a-Mn(a-i)(3-)/2, |5,2(0)| < UWHa-')^ , {3^^ < \\4a-^'+'l (120) 

As was shown above, the operator Catc is self-adjoint in the inner product (|25p and hence {Cabc4>, Cbd) = 
{(j), CabcCbci)- Using (|45p and the fact that QbM is orthogonal to (/>, we obtain the estimate 

\{Cabc4>Xbc^\<\\4{d + af-\ (121) 

Lastly, we estimate (A/", Cbd) which can be written, using integration by parts, in the form 

1 f°° 

{■N'Xbci) = -o / 'l^'^9y{'-fabcy^Cbci)dy, 

^ Jo 
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where, recall, 'jabc is the gauge function (see ([Ml))- Here we used that y<j) -^ as y —^ oo so that the boundary 
terms vanish. Using that 

dyijabcV'^Cbci) = dy{ — -^ e-'^y (bci) = [(— 2 ^ ~ ^ °^)Cb" + y^^dyCb„]jabcy^, 

we find 

\{J^Xbc^)\<\\ic + yT'^4^ + 4ic + yT'^<t>\\'■ (122) 

Estimates p2T|) and (fT22)) give 

1^.(0)1 < (d + a)'"'||</'|| + ||(c + 2;2)-2t^0||2. (123) 

Since we assumed d < aln(a^^), the above estimates imply the following inequalities for fi and r^ 

|/,|<a»-2(lni)*||0||, (124) 

\n\<a'-'\Rom + {aHa-')y\RM\<a'-'[id + a)\\<l)\\ + \\{c + y^r'4'] 
+ iaHa-')Y[\\4 + \\ic + y')-'^4^] 
< {a\n{a-^)Y\\(j)\\ + a'~^\\(j)f , (125) 

The estimates (|120p show that g — -^ — ^(1 + o(l)), provided 

a-^\S„2\, |5oi|, |S'o2|, 1^121, aln(a-i)|5ii|, a|5o2||^ii|, a\Soi\\Si2\ ^ a-Hn{a-^), (126) 

which holds, provided |j(/)|| <C 1. Therefore g is invertible and its inverse is of the form g^^ = infa-M i^^'^i^))- 
Hence the equations (|119p can be rewritten as (|49| - ([50|. with Si{(l),a,b,c) — — and TZi{(j),a,b,c) = — . 
Then the estimates of fi, r^ and g given above, imply (ISTt - ([5^ . 

Proof of estimates ([TTUl) - ([n^ . Use e"''^'/^ = -■J-9ye"°^''/^ and integration by parts to obtain 

1 /-OO 2 

(l,aco)-Y^y Z/(c + y2)2<ie-Vdy (127) 

^ (c2'^ + 4d/ (c + y2)2d-iye-Vdy). (128) 

16a Jo 

To extract the leading part in the last integral above we rescale y — > \/o,y to obtain 
(c + y2)2<i-iye-'^*d2/ - a-^^^ / (ac + y^'^^'ye-'^dy . 

'0 JQ 

Next, split the integral up into 

2 2 /■! 2 /-OC 2 

(ac + 2/2)2''-i2/e-'rd2/= / (ac + 2/^)2^-1 ^^-^dy + / {ac + y^'^-^ye'^^dy . 
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The second term on the left hand side is uniformly bounded in a, d small, so it suffices to investigate the 
first term. Write 

\ac + y2)2d-iyg-^dy= / {ac + y^'^-^ydy + f {ac + y^)^'^-^y{e-'^ - l)dy , 



where again the second term is uniformly bounded in a, d small. Explicit integration in the first term yields 

{ac + y2)2'i-iydy = i-((f + acf' - {acf") . 
4d 

By assumption, there is an e > Oc such that d{a) < a^. In particular, a^ — )■ I as a — >■ so 

(I + ac)^'' - {ac)^'^ = (1 + 2dO{ac)) - (I + 2dO(ln ac)) = 0{d\n -) , 

a 

yielding 

{l,Cbco) = -^il + Oid\n-)). 
Iba a 

The remaining terms are estimated similarly. D 
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